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ABSTRACT 

High-mass X-ray binaries (HMXBs) consist of a massive donor star and a compact object. While several of those systems have been 
well studied in X-rays, little is known for most of the donor stars as they are often heavily obscured in the optical and ultraviolet 
regime. There is an opportunity to observe them at infrared wavelengths, however. The goal of this study is to obtain the stellar and 
wind parameters of the donor star in the X1908+075 high-mass X-ray binary system with a stellar atmosphere model to check whether 
previous studies from X-ray observations and spectral morphology lead to a sufficient description of the donor star. We obtained H- 
and K-Band spectra of X1908+075 and analysed them with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. For the first 
time, we calculated a stellar atmosphere model for the donor star, whose main parameters are: M spec = 15 ± 6 M a , T, = 23/5,kK. 
logoff = 3.0 ± 0.2 and logL/L 0 = 4.81 ± 0.25. The obtained parameters point towards an early B-type (B0-B3) star, probably in a 
supergiant phase. Moreover we determined a more accurate distance to the system of 4.85 ± 0.50 kpc than the previously reported 
value. 
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1. Introduction 


High-mass X-ray binaries (HMXBs) are composed of a com¬ 
pact star orbiting a donor star from which there is an accretion 
of material (see IChatvl 1201 3l . for a review). There is a broad 
literature of the X-ray studies of these types of systems com¬ 
pared to the few detailed studies performed about the donor 
stars using stellar atmosph ere models (e.g. IClark et aD 120021 : 
Gonzalez-Galan et alJl2014l) . The physical characteristics of the 
donor stars harbouring high-mass X-ray binaries, however, are 
essential for a complete picture of the physical processes occur¬ 
ring in these binary systems. 

The binary system X1908+075, also known as_4U 1909+07, 
was first discovered in 1978 bv iForman et alJ ( 197 81 with the 
Uhuru satellite. The source has been observed in surveys car¬ 
ried out with OS07, Ariel, HEAO-1, EXOSAT, and INTEGRAL 
satellites, among others. It is a persistent X-ray source that shows 
flu ctuatio ns of 10% in the soft X-rays, 2-12 keV energy range 
dLevine et al.ll2004h . 


The binary system is a highly absorbed and faint pulsar that 
shows strong X-ray pulsations at a period of 605 s (iLevine et al.l 
120041) . The location of X1908+075 in the pulse-orbital period 
plane (ICorbetl 11986h clearly indicates that this is a wind-fed, 
high-mass X-ray binary. 


Send offprint requests to: S. Martmez-Nunez 
* Based on observations made with the William Herschel Telescope 
operated on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astroffsica de Canarias. 


iLevine et al.l (l2004lf determined the orbital parameters of the 
system (see TableQ} with RXTE/PCA pointed observations. This 
study concludes that the system contains a highly magnetized 
neutron star orbiting in the wind of a massive companion star 
(mass function of (6.07 ± O.35)M 0 ) with an orbital period of 
4.400 ± 0.001 days. Thus, they identify X1908+075 as a high- 
mass X-ray binary. Moreover, this study found changes in the 
optical depth along the line of sight to the compact object corre¬ 
lated with the orbital phase. These changes could be causing the 
observed soft X-ray modulation. 

iLevine et al.l (120041 ) estimate a mass of the donor star in the 
range of 9-31 M e and an upper limit on its size of about 22 R 0 , 
using the mass function and the orbital inclination angle. Pre¬ 
vious values were derived from modelling the orbital phase- 
dependence of the X-ray flux. They also infer a wind mass-loss 
rate for the donor star of > 4x 10 6 M 0 yr -1 , which is larger than 
the expected theoretical value according to IVink et all (20001) . 
Given the high rate combined with the estimated mass and ra¬ 
dius, ILevine et al.l (2 0041 ) concluded that the system might be a 
Wolf-Rayet star with a neutron star companion that could evolve 
and become a black hole-neutron star system in 10 4 to 10 5 yrs. 

[Morel & Grosdidien (2005) analysed near-infrared observa¬ 
tions of stars in, or close to, the error box of HEAO-1/A3. They 
suggest that the optical counterpart of the system might be a late 
O-type supergiant at a distance of 7 + 3 kpc. 

In this paper, we present intermediate-resolution (R ~ 2500) 
infrared spectroscopic observations of the counterpart of the X- 
ray pulsar X1908+075 in the error box reported by the Chandra 
position. For the first time, a non-LTE stellar model atmosphere 
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Table 1. Relevant orbital parameters derived bv lLevine et ali d2004 ) 


Parameter 

Orbital second-epoch Values 

a x sin i [cm] 

(1.43 +0.03) x 10 12 

+90 [MJD] 

52631.383 + 0.013 

^spin ([S]) 

604.684+0.001 

P spin [$ S ] 

(1.22 +0.09) x 10~ 8 

e 

0.021 +0.039 

P orb [days] 

4.4007 + 0.0009 

f(M)[M 0 ] 

6.07 ± 0.35 


analysis ofX1908+075 is performed providing important physi¬ 
cal parameters, such as the stellar temperature (77), the terminal 
velocity of the wind (uoo) and the effective surface gravity (g e ff), 
for the donor star. We also calculate a more accurate value of the 
distance to the binary system and enclose the value of its incli¬ 
nation. We conclude that the optical counterpart is, in contrast to 
previous assumptions, an early B-type (B0-B3) star, probably in 
a supergitant phase. 

In Sect. [2] we describe the details of the observations, fol¬ 
lowed by a short overview of the stellar atmosphere models in 
Sect.0 The results are shown in Sect. [4] discussed in Sect.0 and 
conclusions are given in Sect. [6] 


2. Observations 

In July 2009 (proposal id: 135-WHT45/09A), we obtained 
two intermediate-resolution spectra (K and H-Band spectra at 
55018.97 MJD and 55019.92 MJD respectively) of the opti¬ 
cal counterpart of X1908+075 using the Long-slit Intermedi¬ 
ate Resolution Infrared Spectrograph (LIRIS) mounted on the 
4.2 m William Herschel Telescope (WHT), at the Observatorio 
del Roque de los Muchachos (La Palma, Spain). The instrument 
is equipped with a 1024 x 1024 pixel HAWAII detector. We took 
advantage of the excellent seeing and made use of the 0.65" slit 
in combination with the intermediate-resolution K and H pseu- 
dogrisms. The K configuration covers the 2055-2415 nm range, 
giving a minimum resolving power R ~ 2500 at 2055 nm and 
a slightly higher resolving power at longer wavelengths. The H 
configuration covers the 1520-1783 nm range, giving a mini¬ 
mum R ~ 2500 at 1520 nm. 

The position of the X-ray source was accurately determined 
using the High Energy Transmission Gratings (HETG) on board 
Chandra , as the intersection between grating orders m - ± 1 and 
the zero order image. The coordinates of the X-ray source are 
a = 19 h 10'48.2" and 5 = +07 de s35'51.8", with an estimated un¬ 
certainty of 0.65". We then searched for appropriate candidates 
within the error circle using the UKIDSS database. In Fig. [T] we 
present the K-Band image with the Chandra error circle overlay. 
Only one suitable candidate is compatible with the error circle, 
identified as 2 MASS J19104821+073551 6, in agreement with 
candidate A in Morel & Grosdidierl (120051) . 

Data reduction was carried out using dedicated software de¬ 
veloped by the LIRIS science group, which is implemented 
within IRAffl We use the G2V star HD 182081 and the A0V star 
BD-013649B to remove atmospheric features using the IRAF 
standard procedure. The normalized spectra are shown in Fig. [2] 


1 IRAF is distributed by the National Optical Astronomy Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation 



Fig. 1. UKIDSS image of the X1908+075 region. North is up and east 
is left. The red circle indicates the Chandra position of the X-ray source. 
Size of the image: l'x 1'. 


3. Stellar atmosphere models 


During the analysis, the observed spectra are compared with syn¬ 
thetic spectra calculated with the Potsdam model atmosphere 
code (PoWR). The PoWR code provides a model for a spheri¬ 
cal symmetric star with an expanding atmosphere by iteratively 
solving the radiative transfer equation and the statistical equi¬ 
librium equations in non-LTE. The code further includes energy 
conservation and treats both the wind and photosphere in a con¬ 
sistent scheme. 

The main aspec ts of the code are summariz ed in 
iGrafener et al.1 (120021) and lHamann & Grafeneil (120031) . This 
code has been applied to all kinds of stars that could be potential 
donors in wind-fed HMXB systems, suc h as O and B stars (e.g. 
lOskinova et al.| 201l|: Evans et al.| 2 0121) and Wolf-Rayet (WR) 
stars (e.g. Hamann et al. 20061 : Sander et ahll2012l) . Wind inho¬ 
mogeneities are considered in a so-called “micro-clumping” ap¬ 
proach, assuming that the wind is not smooth, but instead con¬ 
sisting of optical thin c ells with an increased density and a void 
interclump medium (cf. lHamann & Koesterkelil998» . 

The “stellar radius” R „ marks the lower boundary of the 
model atmosphere and is set at a Rosseland optical depth of 
r = 20, where we assume that the hydrostatic equation is ful¬ 
filled. In the quasi-hydrostatic regime we use this equation to cal¬ 
culate a hydrodynamically-consistent velocity field dSander et al.l 
120151) . In the outer part we assume a so-called “beta-law’,’ 


v(r) = Coo(l - y) , (1) 

with Uoo being the terminal velocity of the wind from the donor 
star. Both parts are connected such that a smooth transition is 
guaranteed. That means not only the velocities but their gradi¬ 
ents are equal at the connection point, which is usually close to 
the sonic point where the wind velocity is equal to the speed 
of sound. In the final model for the donor star of XI908+075, 
we adopted an exponent of /? = 1.2, after testing models with 
different (5 values. Models with lower values of /3 down to 0.6 
lead to slightly lower fit qualities. The precise choice of beta 
has hardly any effect on the obtained parameters, for instance. 
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Fig. 2. H and K-Band spectra of the X1908+075 donor star. 
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for the estimated mass-loss rate: reducing [5 by 0.2 requires ap¬ 
proximately an increase of 0.05 dex in M to still reproduce the 
observed strength of the Bry line. 

With a given luminosity L, the stellar temperature T, can be 
obtained from L and R * via the Stefan-Boltzmann relation. Apart 
from L, T, , and v m , the other main model parameters are the ef¬ 
fective surface gravity g e ff and the wind mass-loss rate M. The 
loggeft corrects for full radiation pressure a md , i.e. it considers 
not only Thomson scattering, but also the non-negligible contin¬ 
uum and line accelerations. While the actual calculation for the 
model stratification is depth dependent, we use a mean value of 
Trad := a rad/# to define an output value for g e ff that is related to 
the pure gravitational acceleration via 

logoff = log# + log (l-F rad ). (2) 

The details are described in 1 Sander et all (2,015). In addition, 
the clumping factor D describes the maximum density contrast 
reached in the outer parts of the wind. In our models we assume 
that the density contrast starts around the sonic point and in¬ 
creases outwards until the specified value of D is reached. For 
our O- and B-starjnodels, we assume a solar composition ac¬ 
cording to lAsplund et al. ( 2009 ). see Table[2] An overview of the 
model atoms used in the PoWR calculations is given in Table [3] 

4. Results 

We simultaneously reproduce the normalized line spectrum and 
the spectral energy distribution (SED) of X1908+075 by our cal¬ 
culated PoWR model (see Fig. [3). The SED is compared to avail¬ 
able photometry from IRAC (Sn itzer Science! 2009), 2MASS 
(ISkrutskie et al.fl2006l) and lLasker et all (12008 ). First, the in¬ 
terstellar reddening colour excess £b-v 1? determined. Then, 
we apply the reddening law by iFitzpatrickl ( 119991) to calcu¬ 
late £j_ K and compare this value with the estimated extinction 
curve in the direction towards X1908+075. This curve, shown 
in Fig. |3 was_obtained following the procedure outlined in 
ICabrera-Lavers et~al . ( 200 7): red clump giants are selected over 


a near-infrared colour-magnitude diagram. Since this population 
has a well-calibrated luminosity function, we can compare their 
intrinsic and apparent colour and magnitude to derive both their 
distance along the line of sight and the reddening to which they 
are subjected. Once this (d, £j-K)-curve is established, we can 
compare it with the colour excess of our source (Tj_ K = 1.66) to 
estimate a more accurate dist ance of 4.85 + 0.50 kpc co mpared 
to the previous estimation bv lMorel & Grosdidierl d2005l) . 

Since the total luminosity L is critically dependent on the dis¬ 
tance, we can estimate the luminosity, using our SED fit, once 
we have determined the new distance of the system. We ob¬ 
tain a value of 10 4 XI L 3 . This value is used in our atmosphere 
model. As we already need a model for the SED fit, the whole 
process was iterated until we obtained the consistent solution 
described above. The corresponding stellar radius R „ = 16 R 0 
inferred from temperature and luminosity perfectly agre es with 
the radius deduced for a wind-fed system in lLevine et all d2004l) . 

Prior to modelling the spectrum, we compared our data 
with the standard atlases of OB stars in the H and K bands 
dHanson et alJI 1998112005}) to constrain the parameter space of 
the models. The Hen lines at 1.693 /urn and 2.188 fim are weak 
or maybe even absent as it is hard to distinguish them from the 
noise here. Together with the relatively strong hydrogen absorp¬ 
tion li nes in the H-band, this po ints to a spectral type later than 
09.5 dHanson et all 1 19961120051) . The luminosity class is more 
difficult to ascertain. Several indicators point to a high luminos¬ 
ity: The He I series lines are narrow and the H series can be seen, 
at our resolution, up to Br 19. However, the equivalent width of 
the He I line at 1.700 fim (~ 1.9A) is lower but comparable to 
that of the Brackett (11-4) line at 1.6814 /im (2.8 A). This is 
typically seen in supergiants rather than in main-sequence stars 
where EW(Hel) < EW(Br 11-4). All these properties suggest a 
high luminosity class. In addition, Bry in emission is a spec¬ 
troscopic feature only observed in two objects within the sample 
analysed bv lHanson et al.ld2005l) . Both sources belong to the rare 
class of ON/BN supergiants. 

We calculate a set of stellar models starting with typical 
properties of 09.7-B4 stars. An iterative improvement of the pa- 
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Fig. 3. Spectral energy distribution (SED) of the 4U 1907+09 donor star. The PoWR model (red solid line) is reddened such that it fits the observed 
photometry marks (blue boxes). 



Fig. 4. Estimated extinction curve from the direction_of XI908+075, 
obtain ed following the procedure outlined in ICabrera-Lavers et all 
12007) . 


rameters is performed until a sufficient fit of both, the spectral 
lines and the photometry, is achieved. Figure [5] shows the nor¬ 
malized line spectrum in the H and Ks bands together with the 
best-fitting PoWR model. Most spectral features can be repro¬ 
duced, apart from the observed asymmetric shape in some of the 
hydrogen absorption lines in the H band. The corresponding stel¬ 
lar parameters are compiled in Table[2] Because of the limited 
spectral range and the absence of unblended metal lines, there is 
a certain degeneracy in the spectral appearance that introduces 
a quite large error margin. The observation is best reproduced 
by a model with 77 = 23kK and logged = 3.0, which is in 
line with an early B-type donor. To constrain the temperature, 
we used the Bracket series as well as He i at 1.7 pm and He n at 
1.693 pm. The absence of Hen emission at 2.189 pm excludes 
temperatures higher than ~ 30 kK. In the range between 21 and 
28 kK one can always find a model that reproduces Bry and a 
few other lines sufficiently, but most of them fail to predict the 
strength of the very sensitive He 1 2.059 pm line. As the low ter¬ 
minal velocity deduced from Bry already favours the lower part 
of the deduced temperature range, we found that a model with 
77 = 23 kK best fits the observation. 


For logged, the right wings of the bracket series could be 
used to determine an upper limit and excluded values <: 3.3. 
These stronger effective gravity values would lead to broader 
wings than observed, including an unobserved red absorption 
wing for Bry. 

To reproduce the shape of the observed line profiles, the 
emergent spectrum of the PoWR model was convolved with a 
rotational velocity of lOOkm/s. Given the limited spectral qual¬ 
ity this should be seen as just a rough approximation with a rel¬ 
atively large uncertainty of +50km/s. 

Before the donor was detected, iLevine et all (2004) specu¬ 
lated that the donor might be a W R star. With the first infrared 
spectra of the donor obtained by iMorel & Grosdidierl ( 2005 ) 
it became clear, that this is most likely not the case. Instead, 
iMorel & Grosdidierl ( 2005 ) claim that the donor star would be an 
07.5-9.5If supergiant, based on the emission feature next to the 
He i absorption line at 2.11 pm. This feature is actually a com¬ 
plex structure consisting of He i, N m, C m, and O in emission 
lines. The C, N and O components are originated from high ex¬ 
citation states, and since not all of them are covered in our atomic 
database, this peak cannot be reproduced with the PoWR model. 
While one might argue that this line could be an indicator of en¬ 
hanced nitrogen, this should also affect Nm 2.2471/2.2513 pm. 
Since these lines are not detectable in the spectrum, we keep a 
solar chemical composition. 

The Bry line is the only line in the observed spectra formed 
in the stellar wind. Therefore, the determination of the terminal 
wind velocity can only rely on fitting the wind broadening of this 
line. This fitting results in u m = 500 ± lOOkm/s with a relatively 
large error margin. Significantly larger values would over pre¬ 
dict the emission line width or even lead to additional absorption 
features that are not observed. 

The precise B-subtype is hard to determine from the given 
observations. The relationship between the He I line at 1.700 
pm and the Brackett (11-4) line at 1.6814 pm (the former be¬ 
ing deeper than the latter) clearly shows that the donor should 
be B0.5 or even B0 (Fig. 2 in [Hanson et all 119981) . The ob¬ 
tained temperature of 23 kK, however, favours a slightly later 
subtype around B2. We therefore specify the spectral type of the 
XI908+075 donor star to be in the range between B0 to B3. 
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Fig. 5. The observed spectrum (thin blue line) - H band (left panel), K band (right panel) - is compared with the best-fitting PoWR model (red 
thick dotted line). Noticeable spectral lines are identified. 


Table 2. Derived properties of the X1908+075 donor star. 


Table 3. Overview of the used model atoms. 


Parameter 

Obtained value 

T* 

23.0 + j kK 

T 2 

22.71® kK 

log L/L e " 

4.81+0.25 

En-y b 

5.7 + 0.1 

d c 

4.85 ±0.50 kpc 

Voo 

500+ lOOkm/s 

R, 

16 R q 

log geff 

3.0 ±0.2 

logg 

3.2 ±0.2 

-^spec 

15 + 6 M Q 

M 

lO -655 M 0 /yr 

D d 

10 

^dop 

20 km/s 

Urot 

100 ± 50km/s 

My 

-5.29 

x H e 

0.735 

X H e e 

0.256 

*c e 

2.2 • 10- 3 

x N e 

6.1 • 10~ 4 

Xo" 

5.3 • 10- 3 

We" 

1.1 • 10~ 3 

log Q(/ 

47.1s- 1 


Notes. (a> Luminosity assuming a distance of 4.85 ± 0.50 kpc. ^Ob¬ 
tained via SED fitting, using the reddening law from lFitzpatrickl 11999). 
(cl Distance inferred from consistent reddening fit and extinction curve 
shown in Fig.[4] Maximum value used in a depth-dependent ap¬ 
proach, starting at the sonic point and increa sing outwards. (c) Solar 
abundances, adopted from lAsplund et ahl ( i200 9). m Number of hydro¬ 
gen ionizing photons . 


Ion 

Levels 

Lines 

Hi 

22 

595 

H u 

1 

325 

He i 

35 

0 

Hen 

26 

231 

He hi 

1 

0 

Ci 

2 

0 

Cn 

32 

703 

C hi 

40 

630 

Civ 

21 

703 

Cv 

1 

1 

Ni 

1 

1 

Nn 

38 

496 

N in 

36 

780 

Niv 

38 

210 

Nv 

2 

0 

Oi 

2 

1 

On 

36 

630 

Oiii 

33 

528 

Oiv 

25 

300 

Ov 

2 

1 

Fei° 

1 

0 

Fe n" 

2 

1 

Fe hi" 

8 

25 

Fe iv" 

11 

49 

Fev" 

13 

69 

Fe vi" 

17 

121 

Fe vn" 

11 

52 

Fe viii" 

1 

0 


Notes. (a> A super-level approach is used for this element. Fe actually 
denotes a generic atom that also includes the additional iron group ele¬ 
ments Sc, Ti, MCr, Mn , Co , and Ni. The details and relative abundances 
are given in lGrafener et all 12003) . 


5. Discussion 

5.1. On the accretion scenario 

We are going to demonstrate that the wind mass-loss rate derived 
in our model is compatible with the observed average X-ray lu¬ 
minosity, considering^ Bondi-Hoyle-Lyttleton (BHL) accretion determine the orbital separation of the neutron star to the donor 
scenario dBondi & HovlefT944l and references therein). We esti- star (a). Given o ur spec troscopic mass, the orbital period derived 
mated the fraction of the wind that is gravitationally captured by by iLevine et al.1 d2004l) (see Table Q} and assuming a canonical 
the compact object to estimate the accretion luminosity. First, we neutron star of 1.4M Q , we obtain a = (1.98 ± 0.10) x 10 12 cm. 
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Fig. 6. Relations between the donor mass and neutron star mass. Lines 
of c onstant orbital inclination constructed from the mass function of 
iLevine et ab <2004 1 are labelled. The X-axis range corresponds to the 
spectroscopic mass error margin, while the Y-axis range covers the 
range permitted for a neutron star mass (e.g. lLattimeil2012l) . 


Second, we determine the Bondi-Hoyle accretion radius, and fi¬ 
nally, we estimate the captured fraction of the wind as follows: 


f(a) = 


M captured 
^wind 


r 

accr 

4a 2 


2GMns 

^wind+^orb' 


(3) 


Therefore, the captured fraction of the wind depends on the 
wind velocity and the orbital velocity as a function of the orbital 
separation. From our model we obtain a = 190 + 40 km/s 
for the given a , and we estimate an orbital velocity of v OI b = 
330 ± 20 km/s (considering our spectroscopic mass, a canoni¬ 
cal neutron star of 1.4M 0 and the obtained a), which leads to a 
captured fraction of the wind of f(a) = 0.0044 ± 0.0016. 

Once we know the captured fraction of the wind, we can 
estimate the X-ray luminosity generated by the accretion of this 
fraction of the wind onto the neutron star as follows: 


M/v oo, so an assumption about u m is required to make an es¬ 
timation of M. Since our adopted u m is very low, it is very 
likely that ILevine et al.l ( 2004 ) assumed a higher value of v a ,. 
Therefore, their estimations of M must be corrected by a factor 
F = 500/(;^ cvlne . Even using this correction, it is clear that the 
estimations from the orbital modulation of the Ah yield higher 
values of M than the values obtained in this work. We have used 
a clumping factor D — 10. Decreasing D in order to increase 
M, worsens the quality of the fit. For instance, a more modest 
value of clumping D — 4 (see Appendix [A] in which models 
with different D values are compared), yields a higher mass-loss 
rate M = 4.6 ■ 10 7 ;V/ 0 /yr. In su mmary, there i s a discrepancy 
between the X-rays estimation by ILevine et ali ( 20041) and our 
model estimation (Mp 0 wr)- The plausible reasons are: a density 
enhancement in the vicinity of the neutron star that systemati¬ 
cally increases the measured Ah ; and/or an accumulation of un¬ 
certainties in both methods, including M, D for the stellar atmo¬ 
sphere calculations, and f.T cvlnc , Ah and its variability in X-rays. 

In addition to the comparison with the mass-loss rate estima¬ 
tion from X-rays, we can compare our estimation to the predic¬ 
tions by I Vink et al.l d2000l) for OB stars. These predictions dis¬ 
tinguish two regimes depending on the effective temperature of 
the star (the so-called bi-stability jump). These two regimes are 
also characterized by a different ratio r = Uoc/v esc , where u es c is 
the escape velocity (corrected for the radiative acceleration due 
to Thomson scattering). In the hot regime the effective temper¬ 
ature is over 30kK and r ~ 2.6, whereas in the cool regime the 
effective temperature lies below 22.5 kK and r ~ 1.3. Hence, 
given our parameters (see Table [2}, and v esc = 560 km/s, the 
donor star in X1908+075 belongs to the cool regime. fVink et ali 
( 2000 ) gives an estimation of the mass-loss rate depending on 
the luminosity, the mass, the ratio r, and the effective temper¬ 
ature of the star. Plugging our parameters into the cool regime 
equation, we find Myink — 8 x 10 7 M e /yr. The difference be¬ 
tween our Mp 0 wr and the expected Myink is on the order of V75, 
hinting that a porous wind treatment in the mass-loss predictions 
might possibly solve this discrepancy. 


, _ 6f A7 ns Afcaptured _ GM x sf ( a )Mv,\nd r ( 

f'acci — - g — - g, (4) 

«NS «NS 

where / ; is the efficiency factor of energy conversion and 7 ?ns = 
15km, obtaining a L accr = [(10 ± 3) x 10 36 ] £ erg/s. 

However, we know that the source is persistent in X-rays 
with an X-ray luminosity of ~ (1 - 4) x 10 36 erg/s in the 
2-30 keV band according to lLevine et ali (2004 5 and using our 
distance value. Since we expect significant X-ray emission be¬ 
low 2keV, we estimate the X-ray luminosity in the 0J-50 keV 
energy range assuming the spectral model used by ILevine et al.l 
(120045 . obtaining L x = (1.4 - 5.6) x 10 36 erg/s. Taking an aver¬ 
age value of L x = 3.5 x 10 36 erg/s, we obtain an energy conver¬ 
sion efficiency of 0.36 + 0.13, which is a reasonable value for a 
wind-fed system. 

5.2. On the wind mass-loss rate determination 

Our stellar atmosphere model is used to derive a wind mass- 
loss rate Mp D wR = 2.8 x 10~ 7 M s /yr, significantly lower than the 
value o f > 4x 10~ 6 M 0 yr" 1 derived from X-rays bv lLevine et al.l 
(120041) . ILevine et al.l (12.0041) employed the density stratification 
of the wind expected from a /Maw of the velocity field and 
the continuity equation. Using the orbital parameters that they 
derived, they make a fit of the orbital modulation of the X- 
rays absorption (Ah )■ From this fit it is possible to derive 


5.3. Enclosing the inclination angle of the system 

We can also constrain the inclination of the system using our 
spectroscopic mass (15 + 6 M 0 ) and the mass function given by 
ILevine et al.l ( 2004 ). Figure [6] shows lines of constant orbital in¬ 
clination for the allowed spectroscopic mass range and the per¬ 
mitted neutron star mass range by the equation of nuclear state 
(e.g. Lattimeill2012l) . The allowed inclination values are in the 
range 46° < i < 58°. The canonical neutron star mass of 1.4 M 0 
favours an inclination of ~ 50°. Our determination agrees with 
the estimation of lLevine et al.l (20041) . who used a completely in¬ 
dependent method and obtained the inclination from the orbital 
modulation of N#. 


5.4. On the evolutionary phase 

Our spectrum shows Bry in emission (see Sect. [D, which indi¬ 
cates that the source could belong to the rare class of ON/BN su¬ 
pergiant systems. There are several observational hints that point 
to this classification: 


- This class_ of objects commonly s 
(e.g. |Walborn| T976|) . Comparing our 


show variability 
data set with 
iMorel & Grosdidiei ( 20051) reveals that Bry for XI908+075 
is not only in emission, but varies significantly between 
the two observations. We measured an equivalent width for 


Article number, page 6 ofH 

















































































S. Martlnez-Nunez et al.: The donor star of the X-ray pulsar X1908+075 


VkK 

40 20 10 5 



Fig. 7. Evolutionary tracks using the Geneva Stellar Models 
dEkstrom et al.1 12012 1 for solar metallicity and a rotation rate of 
iijni/iicrit = 0.4 on ZAMS (it corresponds to 170 km/s for an initial 
mass of 20M q ). Initial masses are: 9, 15, 20, and 25M e . The position of 
XI908+075 is marked by the red star. 

this line to be -4.3 A (55018.97 MJD), which differs by a 
factor of 3 compared to the va lue of -1.4 A (52418.54 MJD) 
found bv lMorel & Grosdidierld2005ll . The peak height in the 
normalized spectrum is also about three times higher th an in 
the observation discussed in lMorel & Grosdidie r (2005;) . 

- ON/BN_supergiant systems seem to be nitrogen enriched 
(e.g. lCorti et al.l200~9l) . X1908+075 has likely gone through a 
common envelope phase in the past because of the proximity 
of the neutron star and its companion (R OI b ~ 1.77?,). Conse¬ 
quently, the source could be nitrogen enriched (e.g. lLangerl 
1201 2l) and thu s, show spect r al fea tures similar to ON/BN 
stars as seen in lHanson et al.. ( 2005 ). However, our spectrum 
does not show unblended metal lines, which makes it hard to 
determine the CNO abundances. Models with increased ni¬ 
trogen did not improve the fit and thus we cannot confirm if 
the star is really nitrogen enriched or not. 

- ON/BN su pergiant system s are candidates for interacting bi¬ 
naries (e.g. iWalbornll976l) . The projected rotational velocity 
derived from our datajs high compared to that of isolated OB 
stars dSimon-Dfaz & Herreroll2014l) . This could be the result 
of a tight binary interaction of the system in the past. 

- Figure [7] shows the evolution ary tra cks using the Geneva 
Stellar Models dEkstrom et alj 120121) computed for single¬ 
star evolution and initial masses of 9, 15, 20, and 25 M e . 
Thus, assuming that binary interactions in our system have 
been negligible in the past, we can compare the current posi¬ 
tion of our donor star to these tracks. Indeed, the position of 
X1908+075 in the HR diagram is very close to the 15 Mo- 
track in coincidence with our spectroscopic mass. However, 
as previously discussed, there are several indications that 
point towards a former interaction that would jnake the star 
to appear much younger than an isolated star (iLangerl2012l) . 
Nevertheless, we can neither confirm nor rule out this ex¬ 
pected rejuvenation of the system due to the large uncertain¬ 
ties in our spectroscopic mass and effective temperature. 

5.5. On the asymmetry of some of the hydrogen absorption 
lines 

The hydrogen absorption lines of our data show an additional 
component with respect to the model (see Fig, [5}. These addi¬ 
tional absorptions have a measured blueshift of about 200 km/s 
compared to the rest of the spectrum. The additional absorption 


is comparable to the intensity of the lines in the model. This 
effect is only seen in the hydrogen absorption lines, while the 
helium absorption lines and Bry are perfectly reproduced by our 
model atmosphere. 

A first hypothesis would be that these lines might be shifted 
because of the t idal disruption expected in this kind of close 
binary systems dAbubekerov et al.l 12004tf in which the optical 
companion is pear-shaped, leading to a non-uniform tempera¬ 
ture distribution in the atmosphere of the star. Even though we 
cannot completely rule out this hypothesis since PoWR assumes 
a spherically symmetric wind, it is certainly quite unlikely since 
there is no reason why this effect is only seen in the hydrogen 
absorption lines. 

Another hypothesis is that these additional absorptions are 
due to the presence of a density enhancement caused by the 
interaction with the neu tron star, as seen in the simulations of 
iBlondin et al.l (1199 ll) and lManousakis et al.l (120121) . This hypoth¬ 
esis is supported by the fact that these blueward shifts of the ab¬ 
sorption features are only seen in the H-band spectrum, while 
the lines in the K-band spectrum are perfectly reproduced. Both 
spectra were taken at different orbital phases. The H-band spec¬ 
trum was taken at orbital phase </> = 0.158, with the neutron 
star travelling to quadrature where the projected orbital veloc¬ 
ity towards the observer is maximum. For an orbital speed of 
i; orh = 330 km/s and an orbital inclination of the order of i ~ 50° 
this component would be u OT b sin i = 252 km/s, of the same or¬ 
der than the measured blueshift. Moreover, the K-band spectrum 
was taken during orbital phase </> = 0.374 with the neutron star 
approaching inferior conjunction where the projected orbital ve¬ 
locity towards the observer is zero. As a consequence, the K- 
band spectrum does not show this effect. 

To reproduce the additional absorptions, however, we need 
the column density of the extra absorber to be at least 
equal to the column density from the photosphere to the ob¬ 
server, which is — 1 gcm~ 2 in our model, corresponding to 
Ah ~ 60 x 10 22 at cm' 2 . Moreover, this extra absorber has 
to be hot enough to produce a sign i ficant amount of hydro¬ 
gen lines at level 4. [Manousakis et al.l (1201 2l) simulations corre¬ 
sponding to our mass-loss rate and terminal velocity has shown 
a very modest density enhancement in the vicinity of the neu¬ 
tron star in comparison to our required estimations. Furthermore, 
Levine et al.l d2004l) were able to explain the observed Ah orbital 
modulation appealing to the wind without any further structure. 
Therefore, it is quite unlikely that the observed extra absorption 
is due to the presence of a density enhancement produced by the 
interaction of the neutron star with the stellar wind. 


6. Conclusions 

For the first time, we obtained a reliable set of stellar and wind 
parameters for the donor star in the high-mass X-ray binary 
XI908+075, allowing us to determine various parameters of the 
system. Most of our results are in line with what is expected 
from a wind-fed binary system. The main conclusions of this 
work are: 

- The donor star of X1908+075 is identified with an early B- 
type star (B0-B3). Its main parameters are: M spec = 15+6 M 0 , 
T * = 23*®kK, logg e ff = 3.0+0.2,andlogL/Lo =4.81+0.25. 

- The low effective gravity perfectly fits with a supergiant star, 
favouring this classification. Moreover, the presence of Bry 
in emission could indicate that the source belongs to the rare 
class of nitrogen-enriched B-type supergiant stars. 
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- The source shows a variable stellar wind, since the equivalent 
width of Bry emission line in our K-spectrum is three tim es 
higher than the value found bv lMorel & Grosdidieil(l2005l) . 

- The derived parameters along this work are consistent with 
a wind-fed high-mass X-ray binary within the Bondi-Hoyle- 
Lyttleton (BHL) accretion scenario. 

- The distance to the system is constrained to the value of 4.85 
± 0.50 kpc, which is well in side the previously estimat ed 
range of 7 + 3 kpc reported bv lMorel & Grosdidierl d2005l) . 

- The wind mass-loss rate we have obtained from the spec- 
tral fit is one order of magnitude lower than the lLevine et al.l 
( 2004 ) estimation using X-ray observations. 

- The inclination of the system is constrained, favouring an 
inclination ~ 50°, assuming a canonical 1.4 M e neutron star. 

As a general conclusion, we would like to emphasize that 
to understand wind-fed binary systems in a coherent way, one 
needs to analyse the donor stars with stellar atmosphere models 
to constrain their parameters and to understand how the wind of 
the donor star interacts with its compact companion. 
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Appendix A: On the clumping factor and mass-loss rate dependency 
































